A lossy transmission line as a quantum open system in the standard quantum limit 
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We systematically investigate how to quantize a transmission line resonator (TLR) in a mesoscopic 
electrical circuit in the presence of the resistance and the conductance of the dielectric media. 
Developed from the quantum bath based effective Hamiltonian method for single mode harmonic 
oscillator, the approach we presented in this article is a microscopic theory integrating quantum 
fluctuation-dissipation relation. To qualitatively check the condition under which the TLR can 
behave as a quantum object we study the classical-quantum boundary characterized by the standard 
quantum limit. 
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An ideal one-dimensional transmission line resonator 
(TLR) can be described by a classical wave equation and 
thus can be quantized [1] as usual by following a stan- 
dard procedure - the canonical quantization approach [2] . 
Such a quantization formalism is not treated as a serious 
issue in the usual realities. Most recently, however, the 
situation has been changed radically due to the rapid pro- 
gresses in solid state based quantum computing (QC). In 
one of such QC schemes a one-dimensional TLR is used 
to coherently couple one or more Josephson junction (J J) 
qubits [3-6]. In order to create controllable quantum 
entanglements among these JJ qubits, the TLR has to 
work in a quantum manner as a quantum data bus link- 
ing these qubits. Otherwise the TLR can not induce an 
effective inter-qubit interaction. One can imagine this in 
the conventional cavity QED: the classical cavity mode 
in the strong field limit do not induce inter-atom inter- 
actions to form quantum entanglement of the qubits. 

The above arguments show that the validity of the 
TLR based quantum computing strongly depends on 
whether the mode of TLR is truly quantized, that is, 
has some observable quantum effects. In this sense, the 
way to quantize the modes of the TLR and the corre- 
sponding quantization condition become fundamentally 
important for applications of TLR in quantum informa- 
tion processing. In this article we devote to answer this 
question in a more realistic situation taking into consid- 
eration of leakage. Furthermore, we will try to find out 
what characterize the boundary between the quantum 
and classical regime for the lossy TLR by examining the 
so-called standard quantum limit (SQL) [7,8], which, as 
a consequence of the Heisenberg uncertainty principle, is 
usually referred to as the fundamental limit of the preci- 
sion of repeated position measurements. In the modern 
technology based quantum measurement it was recog- 
nized that if one can reach the SQL in the experiments, 
the quantum behavior is observable even for macroscopic 
objects. Most recently LaHaye et. al. [9] described such 
an experiment with the goal to test SQL on a vibrating 
nano-mechanical beam that is about onc-hundrcdth of a 
millimeter. This excites our interest on the similar prob- 



lem about the actual boundary between the classical and 
quantum regime for the TLR. 

To start with we consider the quantization of the 
modes of the one-dimensional lossy TLR. In our model 
the lossy TLR is treated as an open system interacting 
with a bath, a thermal environment. It may be a back- 
ground electromagnetic field interacting with the trans- 
mitted charge in the TLR, or the classical lead connected 
to the TLR, or other damping mechanism. Mathemat- 
ically the lossy TLR can be well depicted by a wave 
equation with leakage and its Fourier components obey 
the typical dissipation equation. In this sense, the idea 
and methods developed in our previous works [10-12] on 
quantum dissipative system can be applied to the present 
discussion. We will show that the interaction with the 
bath leads to an explicit description for the TLR in terms 
of the generalized Caldirora-Kani (CK) effective Hamil- 
tonian. We also examine how the quantum and thermal 
fluctuations of the environment contribute to the uncer- 
tainty of the canonical variables of the lossy TLR. 

To be universal we firstly revisit the derivation of the 
classical wave function for the lossy TLR in high di- 
mensional case. The model is depicted by four lumped 
parameters, the distributed resistance r, the distributed 
conductance g of the dielectric media, the distributed in- 
ductance I and distributed capacitance c per unit length. 
Let V — V{x, t) and I — I(x,i) be the distributions of 
the voltage and the current vector respectively. With 
the conservation of the current, the Kirchhoff's voltage 
law leads to the equations of motion [13] 



vy = -.i-^g. 



dV 
V-l^-gV-c- (1) 



Eliminating the current vector in the above two equations 
we obtain the high dimensional lossy wave equation for 
the voltage [14]. 

1 d^V 



, dV 

V^V = rqV + {re + gl) ^:- + 
at 



dt^ 



(2) 



where v 



is the velocity of propagation or phase 
speed. For an isolated conductor there does not ex- 
ist the current along the norm direction of the surface 



of conductor and one thus have the boundary condi- 
tion Inlboundary = I ' n\boundary = where n is the di- 
rection of the norm of the boundary surface. Combin- 
ing with the equations of motion, Eq.(l) leads to the 
boundary condition n-\/V\ boundary — n- (— rl — l^^ =0 
for the vohage equation. To be specific in one dimen- 
sional case, we have the boundary conditions for the cur- 
rent /(x = 0, t) = I{x = L, i) = of an isolated 1-d TLR 
of length L. The corresponding boundary conditions for 
voltage equation is Vx{x ~ 0,t) ~ Vx{x ~ L,t) = 0. In 
the following discussion we will focus on the quantization 
of the 1-d lossy TLR with such boundary conditions. 

A canonical quantization scheme for an ideal 1-d loss- 
less TLR can be found in Louisell's monograph with an 
implied Hamiltonian [1]. But for the lossy TLR as an 
open system the energy does not conserve and thus there 
does not exist a time-independent Hamiltonian a priori. 
Hence it is necessary to develope a quantization scheme 
only based on the classical equation of motion. In the 
very original paper of quantum mechanics by Hcisen- 
berg, however, without knowing the Hamiltonian or La- 
grangian beforehand the classical equation of motion is 
sufficient to solve the quantized energy levels of linear os- 
cillator. The modern version of this idea was presented 
in the famous textbook by Landau and Lifshitz [15]. Fol- 
lowing this method we even quantized the quantum dis- 
sipation system [10-12] to re-deduce the so-called CK 
effective Hamiltonian [16-18] directly from the classical 
equation of motion for the open system with the linear 
coupling to the Ohmic bath [19]. 

Now we quantize the modes of the lossy TLR by devel- 
oping our previous approach mentioned above. We can 
start from the lossy voltage equation 

il9 + rc^i-V + rgV-^=0 (3) 



P„(t) = M„e^*-K(t) 



(6) 
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:V + (Iq + rc:^ —V + rgV 
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v^ dt^ dt ' dx^ 

with boundary condition Vx (0,t) — Vx {L,t) — 0. The 
time-dependent Fourier's components Vn (t) of V{x, t) 

corresponding the normal mode i/-|-cos^^ just satis- 
fies the dissipation equation 

^K(t)+7|KW+^XW = (4) 

where the dissipation rate and the effective frequency are 
7 — v^{gl + re) and oj'^ — v'^[rg + (■^) ] respectively. 

To find out the canonical commutation relation for the 
dynamical variables Vn (i) and -^Vn{t) we calculate the 
time evolution of their commutator B{t) = [Vn, -^Vn (i)] 
using the equation of motion (4), obtaining the close 
equation -^B = —^B. That means B oc e~^*, or 

[K, jVn (i)] - ihe-^'/Mn. (5) 

Here we choose the effective mass Af„ = -%- so that we 
can get a correct expression for the energy of TLR. There- 
fore, we can define the canonical momentum operators 



to realize the conventional canonical commutation rela- 
tion [Vn,Prn] = «M„„. 

In terms of the annihilation and creation operators 
a„ [t] and ajj (t) defined by 



Vn{t) 
Pn{t) 



-^ {an + al_ 



(7) 



the C-K Hamiltonian for the 1-D lossy TLR is obtained 
as 7Y = Y^n^n, where 

Hn = fiLOn{[a!^ sinh(7t) + al^an cosh(7i)] + h.c}. (8) 

Here we have ignored a time-dependent c-number 
hi^n cosh(7t) 

The above effective Hamiltonian Ti can force the TLR 
to evolve into the multi-mode squeezed state when the 
TLR is initially prepared in a multi-mode coherent state 
\a) = |ai, a2,..., a„, ...) where \a)^ denotes the n-th 
mode. This conclusion can be proved by rewriting i/„ 
as 



H = nuj[^AlAn 



Sl (t) haJna^anSn (t) 



(9) 



Here, {An} is a new set of bosonic operator defined as 
the unitary transformation of a„ An = Sj^ (t) a„5'„ (t) by 
the squeezing operators [20] 



Sn (t) = exp 



^7^(4^ - al) 



(10) 



We can also write down the explicit expressions An = 
^cin — ?7ajj with ^ = cosh(7i/2), 77 = — sinh(7t/2). Then 
we obtain the evolution operator 



Un {t) ^ Sl (t) cxp[-iu!nalant]Sn (0) . 



(11) 



which gives the the final state |^„ (t)) = [/„ (t) |a)„ — 
|Q,g-»w„t^ ^^ _j^^ g^g g^ squeezed state defined as the eigen- 

statc of An with eigenvalue ae"'""*. 

Studying the above results carefully, it seems that the 
above results are not totally convincing due to the viola- 
tion of the uncertainty principle about coordinate Vn (t) 
and momentum P„ (t) because the above arguments are 
too phenomenological and the source of dissipation is 
not considered microscopically. Thus the Brownian mo- 
tion can not be analyzed in the frame of the effective 
CK Hamiltonian formalism. What's more, if we would 
use the above CK Hamiltonian without restriction, some 
ridiculous conclusions are to be reached. Fortunately, we 
can solve this problem by demonstrating the derivation 
of the phenomenological CK Hamiltonian, starting with 



the conventional system-plus-reservoir approach. Sup- 
pose the environment is a bath of many harmonic os- 
cillators linearly coupled to the open system. The bath 
and system constitute a conservative composite system 
and thus its quantization is rather straightforward. As 
shown in our previous works, the total wave function 
is partially factorized with respect to system and bath 
when the Brownian fluctuation can be ignored under cer- 
tain conditions [10]. the factorized part of the system is 
just the CK effective wave function governed by the CK 
Hamiltonian. With the these recalls we now consider the 
damping mode equation 



Vn (t) +lVn (t) + UjlVnit) - /„ (t) 



(12) 



where the parameter 7 and a;„ are the same as before 
and 

]n\t) = - > Cj„(a;„jo costJ„ji + x^^^ ) (13) 



is a Brownian driving force. Here, x^o (ijo) are the ini- 
tial values of the canonical coordinate (its first deriva- 
tive with respect to i) of the harmonic oscillators of the 
bath that coupled to the transmission line, LOnj is the fre- 
quency of the j-th oscillator of the independent reservoir 
coupling to the n-th mode of QTL with the correspond- 
ing coupling constant Cj„. 

The solution Vn {t) = Qn (t) + X^i ^nj (t) of the above 
motion equation (12) can be solved as a direct sum of the 
dissipativc motion 

Qn (t) = ai (t) VnO + 02 {t) T4o (14) 

and the quantum fluctuation of the reservoir. 

£,n] (t) = 2_, ^"Jl (^) ^"i" + ^"J2 (*) ^njO (15) 



where the coefficients ai (t) 



^(i^e-^* - fie-"*) 



a2{t) = T;^{e ^*-e •"*) with fi = ^ +iuj„, i^ = ^ -luj. 



X 

2 



— 2 

2 ""ni 



\ — ^ and we ignore the lengthy and too te- 
dious expressions of bnji{t) and 6„j2 (t) ■ 

To revisit the quantum effect of QTL, we need to con- 
sider the standard quantum limit (SQL). If the measure- 
ment that probes it can reach the accuracy of the SQL, 
then the quantum effect is observable. Obviously the 
standard quantum limit of Vn is contributed by the mea- 
surement of QTL and the bath fluctuation 

(AV;««^)2 = {AQnf + al (16) 

where AQn is the variation of Qn(t) and 

al=J2{A^n,r (17) 



is just the width of Brownian motion. The first part of 
the above equation is just an average over the pure quan- 
tum state of the TLR. Thus we can use the uncertainty 
principle to determine the value of SQL. According to 
the commutation relation Eq.(5), the uncertainty rela- 
tion means 



AK > 



?iexp(— 7t) 

2M„AK 



(18) 



Then 



(AQ„)' > 



2h 



MnUo'^ 



■e"27t [^ sin2 J^t + oj'n sin Sw^i] (19) 



It is easy to see that as time t ^ cx), the standard quan- 
tum limit jAQ (t)| ^0 (see Fig.l). If the quantum fluc- 
tuation caused by the bath fluctuation ^„j (t) is ignored 
inappropriately, the standard quantum limit is zero af- 
ter a long time evolution! We will show as follows that 
the quantum fluctuation contributes a nonzero part as 
compensation. 




FIG. 1. The time evolution of the standard quantum limit 
of the system in low temperature limit. The solid line shows 
the evolution of the SQL of the TLR while the dashed line 
shows that of the environment. Here the time is scaled in the 
unit of - 

7 

We notice that the second part of Eq. (16) is the ther- 
mal average over the reservoir states and the fluctuation 
On at temperature T is 



ai(t)^ycoth(^ 

"^ ' '^ \2kBT 



n[h%At)+^ljbl32{t)\ 



ZTHfij UJfij 



(20) 



where ks is the Boltzman constant and T is the temper- 
ature. This is just the width of the Brownian motion, 
which characterizes the extent of fluctuation around the 
damping path Qn- In the classical limit, known as the 
"Ohmic friction" condition, the fluctuation can be eval- 
uated as 



^lit) 



hj 



TlMnUj'n JO 



dcj coth 

ft ^ sin wit 



(1 -2e"-(7 
—2uj sin iu'^t sin cot) 



\2kBTj "■ ' 

+ u'n COS Uj'nt) COS LOt 



(21) 



where L{lo) = w/[(w^ - lu"^) + "/"^lu"^]- 

As can be seen from Fig.l, in the low temperature 
hmit cr^(i) is zero initially, and then approaches its final 
equilibrium value in a time interval of the order of I/7. 
The above results shows the important limit on the stan- 
dard quantum limit of TLR caused by the environment 
fluctuation. It seemingly depends on the details of the 
reservoir, but they can be universally summed up to the 
observable quantities of the TLR in certain limit case. 
When the time is large enough, the total fluctuation of 
the whole system approaches (see Fig. 2) 



(AV„ 



'SQL 



it 



2ttM^ujL 



arctan 



0J„ 



1< 



(22) 




FIG. 2. The time evolution of the standard quantum limit 
of the whole system at low temperature limit. Again the time 
is scaled in the unit of - 

If the damping rate 7 is much smaller than the fre- 
quency of the oscillator a;„, this width happens to be- 
come the same as the width h/{2MniOn) of the ground 
state of the mode, which is just the SQL of the quadra- 
ture amplitudes of the oscillator mode Vn (t) . 

To conclude this article let us estimate the above re- 
sult numerically according to the parameters given in the 
experimental proposal in ref [3,4]. The eigenfrequency of 
the TLR mode in resonant with the Josephson junction 
qubit is about 10 GHz while the dissipation rate is about 
6.25 X 10*^ Hz. Then 



a(k: 



SQL 



0.2fiV 



where L is the length of the TLR. Thus if the precision of 
the experiment can reach this limit, the quantum effect 
can be observed. We also note that our exploration in 
this article reveals the close relation between the SQL of 
the open system and the quantum fluctuation of the en- 
vironment: starting from the initial semi-classical state, 
each mode of TLR experiences a damping squeezing (see 
Eqs.(9-ll))bcfore it reach the SQL. Once near the the 
SQL the quantum fluctuation take place to play as a 
quantum noise against the infinite squeezing by quantum 
dissipation. 



This work is supported by the NSFC and the Knowl- 
edge Innovation Program (KIP) of the Chinese Academy 
of Sciences. It is also funded by the National Fundamen- 
tal Research Program of China with No. 001CB309310. 
One of the authors (CPS) thanks Y. B. Dai for useful 
discussions. 



[2] 
[3] 

[4] 

[5] 

[6] 

[7] 
[8] 



[9] 
[10] 

[11] 
[12] 

[13] 



Electronic address: suncp@itp.ac.cn 
Internet web site: http:// www.itp.ac.cn/~suncp 
W.H. Louisell, "Quantum Statistical Properties of Radi- 
ation", John Wiley and Son's, New York, (1990). 
R A. M. Dirac, Proc. Roy. Soc. A133, (1931) 60. 
Alexandre Blais, Ren-Shou Huang, Andreas WallrafI, S. 
M. Girvin, R. J. Schoelkopf, cond-mat/0402216, 2004. 
S. M. Girvin, Ren-Shou Huang, Alexandre Blais, Andreas 
WaUraff, R. J. Schoelkopf, cond-mat/0310670, 2003. 
R. J. Schoelkopf, A. A. Clerk, S. M. Girvin, K. W. 
Lehnert, M. H. Devoret, cond-mat/0210247, 2002. Also 
in "Quantum Noise in mesoscopic physics" (Kluwer Ac. 
Publ., 2003) 

S. Camalet, J. Schriefl, P. Degiovanni, F. Delduc, cond- 
mat/0405597, 2004. 

L. D. Landau, R. Peierls, Z. Phys. 69, 56 (1931). 
V. Braginsky and F. A. Khalili, "Quantum Measure- 
ment", Cambridge, 1992; V. Braginsky and F. A. Khalili, 
Rev. Mod. Phys, 68, 1 (1996). 
M. D. LaHaye et al.. Science, 304, 74 (2004). 
L. H. Yu, C. P. Sun, Phys. Rev. A, 49, 592 (1994). 
C. P. Sun, L. H. Yu, Phys. Rev. A, 51, 1845 (1995). 
C. P. Sun, Y. B. Gao, Phys. Rev. E, 49, 592 (1998) 
To examine it in the two-dimensional case we consider 
the Green's theorem 



I ■ dl = (/) (V • I)ds 



where dl —{—dy,dx) along the direction of tangent. The 
close path integral in the left hand side represents the all 
charges going into the domain surrounded by the closed 
path, which is just the sum of all leakages within this do- 
main, we can consider the integral in the bulk cell Ax Ay, 
leading to the Eq.(l). 
(23) [14] For the current vector, however, instead of a lossy wave 
equationm the modified equation of motion of current 
vector is 



aj 



d^3 



[15] 



[16] 
[17] 



V (V ■ J) = RG3 + (Rc+IG) — + Ic— 

But in one dimension, it can be reduced to the generic 

lossy wave equation. 

L. D. Landau and E. M. Lifshitz, "Quantum Mechanics 

(Non-relativistic Theory)", 3rd Edition, Pergamon Press, 

1987. 

P. Caldirola, Nuovo Cimento, 18, 393,(1941) 

E. Kanai, Prog. Theor. Phys. 3, 440 (1948) 



[18] H. Dekker, Phys. Report, 80, 1 (1981) 



[19] A. O. Caldeira, A. J. Leggett, Ann. Phys. 149, 374 (1983), 

and Physica 121A, 587 (1983) 
[20] Y. D. Wang, Y. B. Gao and C. P. Sun, cond-mat/ 

0212088, to be published in Euro. Phys. Journal B. 



